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Abstract 

We consider Andreev reflection in a two dimensional junction between a nor- 
mal metal and a heavy fermion superconductor in the Fulde-Ferrell (FF) 
type of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state. We assume s- 
wave symmetry of the superconducting gap. The parameters of the super- 
conductor: the gap magnitude, the chemical potential, and the Cooper pair 
center-of-mass momentum Q, are all determined self-consistently within a 
mean-field (BCS) scheme. The Cooper pair momentum Q is chosen as per- 
pendicular to the junction interface. We calculate the junction conductance 
for a series of barrier strengths. In the case of incoming electron with spin 
o =f only for magnetic fields close to the upper critical field H c2 , we obtain 
the so-called Andreev window i.e. the energy interval in which the reflection 
probability is maximal, which in turn is indicated by a peak in the con- 
ductance. The last result differs with other non-self-consistent calculations 
existing in the literature. 
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1. Introduction 

Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) superconducting state has been 
proposed theoretically in the early 1960s [H, 0]. m this unconventional su- 
perconducting state the Zeeman splitting for electrons at the Fermi surface 
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makes it favorable for the Cooper pair to have a nonzero total momentum 
Q = 2q and consequently, the phase of the superconducting gap parameter 
oscillates spatially with the wave vector Q. By forming such condensate of 
moving Cooper pairs, the superconducting state survives to magnetic fields 
higher than the Pauli H C 2 limit. The FFLO state has suddenly gained re- 
newed interest recently because of its possible realization in the heavy fermion 
superconductor CeCoIn 5 jijl, HJ, as well as in the layered organic supercon- 
ductors k-(BEDT-TTF) 2 Cu(NCS) 2 0] and ^'-(ET^SFsCHaCFaSOa 0. All 
those systems have a reduced dimensionality, what is crucial for FFLO phase 
stability, as then the orbital effects are suppressed and the Pauli effect (Zee- 
man splitting) may become the dominating factor. The FFLO state is also 
investigated in context of its possible realization in high density quark and 
nuclear matter [zj], as well as in optical lattices jsj. 

The irrefutable evidence for the FFLO state should be based on the phase 
sensitive experiments, because they can reveal the spatial variation of the 
phase or the sign of the order parameter. One of such experiments is the con- 
ductance spectroscopy in a normal metal (N) - superconductor (S) junction 
(NSJ). The method for analysis of the NSJ conductance has been provided by 
Blonder, Tinkham, and Klapwijk j9j. Such experiments have already turned 
out to be successful for example in the determination of the gap parameter 



d x 2_ y 2 symmetry in the CeCoIn 5 system [10j. A crucial role in the conduc- 
tance spectrum is played by the Andreev reflection (AR) processes (llj . In 
the simplest view of the Andreev reflection an incident electron entering from 
N into S is converted at the NSJ interface into a hole moving in the opposite 
direction and Cooper pair inside SC. 

The conductance characteristics for a NSJ with superconductor in the 
FFLO state has already been investigated for both the cases of FF (with 
A(r) = A Q e iQr ) H and LO (A(r) = A Q cos(Qr)) [13] types of the FFLO 
state, as well as for the case of superconductor with supercurrent [13, Hi 



(i.e. the situation similar from formal point of view). The emphasis in these 
papers is put on the case with the d-wave symmetry of the gap (see also 



16l 1 171 Il8l . Il9| for the case of a BCS state with <i-wave symmetry). 



Here we consider an s-wave superconductor in the FFLO state of the 
FF type with the qvector oriented perpendicular to the junction interface 



(similarly to Ref. [13[ and different from Ref. [12| in which q is parallel to 
the interface). Recently, similar investigation has been performed for a one- 
dimensional situation 20| and it has been shown that with superconductor 
in the FF state there are lower and upper bounds on energy of the incom- 
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ing electron for the Andreev reflection to take place (an Andreev window 
appears). Aim of the present paper is to analyze the Andreev window in a 
self- consistent model. Namely, we choose such Cooper pair momentum 2q, 
which minimizes the free energy of the system (contrary to Refs. 13|, |20|], in 
which q is fixed) and we adjust the chemical potential fi so that the particle 
number n is kept constant. We show that such careful examination of the 
superconductor properties leads to important alterations of the conductance 
spectrum. 

We analyze the situation for heavy quasiparticles, because the FFLO state 
is more likely to appear when the carriers are heavy (then the orbital-effect 
influence is suppressed). For heavy fermion systems it is peculiar why AR 
takes place at all. Based on the BTK theory due to a large Fermi-velocity 
mismatch it should be severely limited by a high effective barrier strength Z. 
Nevertheless, AR is observed in those systems and theoretical efforts have 



been made to understand why it is the case [22|, [23|, [24|. Here we neglect 
the Fermi velocity mismatch by assuming equal masses of quasiparticles and 
equal chemical potentials on both sides of the junction. 

The structure of the paper is as follows. In Section [2] we discuss the super- 
conducting state of quasiparticles with heavy effective masses in both two- 
(2D) and three-dimensional (3D) cases. In Section [3] we present the results 
concerning Andreev reflection in a normal metal - heavy fermion supercon- 
ductor junction. Finally, in Section H] we provide a summary and outlook. 

2. Fulde-Ferrell superconducting state characteristics 

We consider here 2D and 3D system of paired quasiparticles with heavy 
masses m = 100 m e , which roughly corresponds to the heaviest band of 



CeCoIn 5 [25(. We assume s-wave pairing symmetry, which has not been 
analyzed in detail as yet, since most of condensed matter correlated-electron 
systems hosting FFLO state are believed to have a <i-wave gap symmetry. 

The system of self-consistent equations describing the superconducting 
state is very similar to the one presented in [26[ (see Section V) except that 
here we do not take into account two characteristic features for strongly- 
correlated electrons: the spin-dependence of quasiparticle mass and the ef- 
fective field acting on them. These features appear in both slave-boson and 
Gutzwiller approaches to correlated electrons, but in the latter only if the 
Gutzwiller band-narrowing factors are taken as dependent on the system 
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spin-polarization and supplemented with Maximum Entropy based treat- 



ment 



27j . In this respect, we analyze here the more standard situation and 
the case with spin-dependent masses and effective field is left to a separate 
study 

Explicitly, the free energy functional, the gap parameter Aq for the 
center-of-mass- momentum Q, and the condition for the number of particles 
n (per atom) have the respective forms: 

T = -^TVlna + e-^+E^-^O+iV^ + ZiiV, (1) 



kcr k 

V ^l- f(E^) - f(E H ) 



Aq - W A Q , (2) 

k k 



n 



n t + H = £ Yl {<f( E ^) + ^[1 " /OV,)]}, ( 3 ) 

kcr 



where J-(T, H a ; n, A q ) is the system free-energy functional for the case of 
a fixed number of particles [28| (we set the band filling n = 0.97), V is 
the interaction potential, u^, are the Bogolyubov coherence coefficients, 
f{E\ ia ) is the Fermi distribution, and n a is the spin-subband filling. The 
physical solution is that with a particular Q which minimizes ([I]). The state 
with Q = is called the BCS state, and that with Q ^ - the FF state. 
The quasiparticle spectrum is characterized by the energies (c.f. also (29| ) 



£ kc r = E k + a&\ £ k = ^k S)2 + A Q , (4) 

£k } = 2^ k+qt + £- k+c ^)' & = 2^ k +qt ~ C-k+q4-)> (5) 
with the dispersion relation taken in the computation as 

h 2 k 2 

6- = "2^- -<rh-», (6) 

where h = gfiBH a and H a is the applied magnetic field. Exemplary phase 
diagrams obtained for these parameters on the applied field H a and temper- 
ature T plane are exhibited in Fig. [TJ 

3. Andreev reflection: self-consistent solution 

For the analysis of the Andreev reflection process we take the parame- 
ters obtained self-consistently (from the procedure presented above) for the 



4 



5 
4.5 

4 
3.5 

3 
2.5 

2 
1.5 

1 

0.5 




Phase Diagram - 2D 



Momentum (1/A) 
0.014 



NS 



BCS 



Phase Diagram - 3D 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1 
Temperature (K) 



Momentum (1/A) 
0.016 




1 1.5 
Temperature (K) 



Figure 1: (Colour online). Applied magnetic field-temperature phase diagram for a 2D 
(a) and 3D (b) system of paired heavy quasiparticles. For further analysis of the Andreev 
reflection we take the parameters obtained along the T = Q.02K s» line. Note the 
extended regime of the FF state stability in the 2D case. This is only due to geometrical 
reasons, as we disregard orbital effects. 



superconducting state. Although this state is either 2D or 3D we consider 
only 2D NSJ for simplicity. Kinematics of the reflection may be analyzed by 
means of the Bogolyubov-de Gennes (BdG) equations 

(i'W -£$ r ))«r,*) = ^(r,,), (7) 
where the two-component wave function is given by 

«r,,) S (8) 

where a = ±1 is the spin quantum number and a = —a. The single-particle 
Hamiltonian is the following 

ffo*(r) = -&h-(i+ U(r), (9) 

where r = (x, y) and the interface scattering potential is chosen as a delta 
function of strength if, i.e. U(r) = H 5{x). 

We choose the Fulde-Ferrell (FF) type of superconducting state, in which 
A(r) = Aqe* 2qr and set the direction of the Cooper pair momentum Q = 2q 
as perpendicular to the junction interface (similarly as in Refs. 13| and 201]). 
namely Q = (Q, 0). The junction geometry is presented in Fig. [2j The plane 
wave ansatz 

Mr,*) = e te ( ^X'^) ) > ( 10 ) 
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Figure 2: Junction geometry for incoming particle of spin a =j. Normal-state and Fulde- 
Ferrell regions are marked. Interface lies at the x = line. Full circles mark quasiparticles 
and empty ones mark quasiholes. Momentum of each of them is marked with a boldface 
letter, and amplitude with an italic letter. Namely, incoming particle has the momentum 
k, and amplitude 1, reflected hole has p and a, reflected quasiparticle: k', b, transmitted 
quasiparticle: kj , c, and transmitted quasihole: kj, d. The angle of incidence is equal 
to 9 and to the angle of reflection but other angles (of reflection of quasihole and those of 
transmissions) may differ by a small amount (c.f. also Fig. 2]). 



with u and v as constants, leads to the following matrix equation 
-£ + 6c +q , CT A Q \ ( u\a) 



Aq -E - 6c-q,a J \v\a) 



0, (11) 



where quasiparticle energies are given by flS]). Equation (ITT]) gives the dis- 
persion relations for quasiparticles and quasiholes in the superconductor 



E k± = { ^ Q' 1 ( 12 ) 

where £^ have been defined earlier in (jSJ). One may check that the above 
equation is in accordance with (TjJ, as E k+ = E^ (quasiparticle) and E k _ = 
— E^i (quasihole) for incoming particle with spin a =y, as well as E k+ = -EL14 
(quasiparticle) and E^_ = — -E-kf (quasihole) for incoming particle with spin 
a =] r . The plots of the energies ( lT2"j) are exhibited in Fig. El 
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Figure 3: (Colour online). Quasiparticle and quasihole dispersion relations E\^±: (a) for 
BCS state at H a = 1.01 T; (b), (c) for FF state at H a = 4.61 T < H c2 (b) and H a = 4.01 T 
(c) for incoming electron with spin a — 1 (b) and er = — 1 (c); (d-f) closer view at the 
region around k sa +fci? with magnetic field close to H c n. (d) shows E\^± for a set of 
magnetic field H a values; (e) and (f) show Ek± for a choice of angle of incidence 9 values 
for 2D (e), and 3D (f) superconducting states. The angle of incidence is = in (a)-(d). 
The quasimomenta kf 2 marked in (a) and (b) are solutions to the equations E = E^± 
propagating in the positive (superscript "+") and negative ("— ") x direction. Situation 
with incident particle with energy E in (a) (E' in (b)) results in Andreev reflection process. 
Note that the dispersion relations for BCS (a) do not depend on the incident angle 8. 
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As we consider electron injected from the conductor side of the junction, 
the corresponding wave functions can be expressed as (we have omitted the 
spin part for clarity) 



<J<(r) =^J) e- k ' + a ( J ) e*" + 6 ( J J e' kV , (13) 

*> W = d ( ) + C ( ; 2 ?<- ) (") 

where i[)<(t) an d V ; >( r ) describe wave function on the normal-metal and 
superconductor sides, respectively. The quasimomenta kf (for quasihole) 
and (for quasiparticle) are solutions of ( !T2|) for a given incident energy 
E propagating in the positive x direction. From the translational symmetry 
of the junction along the y direction comes conservation of the y momentum 
component. Namely, k y — k' y — p y — k\ y = k^y All the wave vectors are 
presented in Fig. 4. 

We use standard boundary conditions 

V'<( r )U=o = '0>(r)| a >=o, ( 15 ) 
dtp < (r) l 9^>(r), 2mH 

x=o Z5-lj><( T )\x=0- ( 16 ) 



O \X = r% \X = V a-O 

ox ox h z 

Those conditions lead to the following set of 4 equations for the amplitudes 

(a, b, c, d) 

1 + b — cu 2 — du\ = 0, (17) 
a — cv 2 — dv\ = 0, (18) 

2777 H 

ik x {l - b) - cu 2 i{q + k+ x ) - dmi(q + kf x ) + + b) = 0, (19) 

nr 

, , 2mH 
aip x — cv^iyk^x — q) — dv\i(ki x — q) H ^— a = 0, (20) 



which are very similar to those in e.g. [30], but for our case vectors are 
replaced by their x-components: e.g. k ■<->• k x , p -<-> p^. From its solution one 
can obtain probabilities of hole reflection p a rh = |a| 2 ^jpi, particle reflection 
Pre = l^l 2 ' quasiparticle transmission 

- _ Lia lW-Wjl + g , 91 x 
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a) 



Normal state 



BCS Superconductor 




b) 



Normal state 



FF Superconductor 



E =E 




Figure 4: (Colour online). The junction geometry in the reciprocal space. All vectors are 
marked. It can be seen that only the incident and reflection angles are equal to 9. It can 
be anticipated at this point that changing 9 for BCS state does not lead to drastic changes 
in the transmission/reflection probabilities, whereas for the FF state the situation is quite 
different since q ^ induces anisotropy in the reciprocal space. The energy E value has 
been chosen as 10 K for all graphs except (b) "FF Superconductor" for which E = 0.01 K 
(for E > 0.5 if there would be no E = i?k- regions in this case). The dashed lines are 
guide to eye and illustrate the conservation of momentum y-componcnt. 



9 



and quasihole transmission 



(t . , JI 2(l«il 2 -KIWi + J + ? 



ps, = ( 22 ) 



where the a superscript indicates the spin of the incoming electron. These 
probabilities are plotted in Fig. |5]for various barrier strengths Z = 2mH/ (kph 2 ), 
where we define Fermi wave vector kp using the zero- field value k,F = \\f2mp,. 
Note also that we do not use the assumption k = k' = p = kf = k£ ~ kp 
utilized at this point in majority of papers on Andreev reflection, because we 
deal with heavy quasiparticles for which \x is of the order of 100 K. There- 
fore the assumption \i 3> E is not, strictly speaking, applicable in the present 
situation. 

As can be seen in Fig. there are lower and upper bounds on the energy 
of incident electron for the Andreev reflection to take place. For example in 
(b) this region is [0.02 K, 0.1K] and corresponds to the region in Figs. [3b 
and EH between the maximum of Ek- and minimum of E^ + . The position and 
existence of the Andreev window is strongly dependent on the parameters 
chosen. Namely, a deviation from H c2 by 1% (for a 2D system) or from 
the perpendicular incidence by 5 (2D) or 28 (3D) degrees makes it disappear 
(the lower bound on the energy E at which the Andreev reflection takes place 
goes below E = 0). It is clear that the Andreev window persists to larger 
degrees for a 3D system (for geometrical and energetical reasons, c.f. also 
Fig. |3f). For incident electron with spin a = — 1 the situation is different as 
Andreev window appears for high E values (e.g. between 5.2 K and 6.1 K 
- see Fig. [3b and Fig. \5jp). Note also the qualitatively different behavior of 
the probabilities for the BCS (a) and FF (b-d) superconductor. Namely, the 
reflection of a hole probability p r h for FF is constant in the whole Andreev 
window region and decreases less rapidly with the increasing Z (compare the 
middle graphs of Fig. [5] (a) and (b)). This behavior produces very different 
conductance peaks for the two cases. 

Differential conductance (G = dlldV) can be obtained from the reflection 
and transmission probabilities 0, [2l| in a straightforward manner 



Gns = l 2 d6cos6[l -p° re (E 1 9)+p; h (E,e)]. (23) 

1 J-tv/2 



The final result of our calculation is the conductance G averaged over spin 
and normalized with respect to the conductance G a nn of the junction with 
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Figure 5: (Colour online). Probabilities of reflection and transmission processes (a) for 
BCS state and several barrier strength Z values; (b) for FF state for barrier strengths 
Z = (top), Z = 0.5 (middle) and magnetic field H a = 4.51 T (bottom); (c) for FF state 
with incoming electron of spin a = — 1 and field H a = 4.01 T; (d) for 3D FF state for 
several angles of incidence 9. The angle of incidence is taken 6 = unless stated otherwise. 
Note that the probabilities for BCS state (a) do not depend on 9. 
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A = 0. Namely, 

G = ™ ™ . (24) 

This quantity is exhibited in Fig. EJ In Fig. Efc also the spin-resolved 
conductance G a = G a ns jG a nn is presented. 

The Andreev window for o =f manifests itself by a peak at E — 0.02fT 
(b) and E = 0.5K (d). This peak is not present in 2D case for H a < 4.535 T 
(see (c)) because then the quasihole energy E^ falls below E = (see Fig. 
E]) and there is no Andreev window for a =|" in this case. Note that in the 
non-self-consistent calculations [lil, [io[ the a =t peak was fully present (c.f. 



Fig. 4a of Ref. [13|). Note also that (a) agrees qualitatively with the results 
of e.g. Ref [12]. As should be expected the a =f peak is more pronounced 
for 3D system, but for 2D it is also present for the applied field close to H c2 . 

From the comparison of (a) with (b-d) it is evident that for FF state the 
conductance peaks are much broader than in the BCS case. This is because 
of the anisotropy in the reciprocal space induced by q ^ (see Fig. @J, which 
gives ^-dependent position of Andreev window. Therefore, integration of the 
[1 —p° e (E, 6) +p° h (E, 6)} factor over 6 in ( 123]) gives broader peaks in the FF 
case. 

4. Conclusions 

In this paper we have provided a detailed analysis of the Andreev reflec- 
tion from an FFLO superconductor (of the FF type) by analyzing it within 
a fully self-consistent scheme for the s-wave symmetry of the superconduct- 
ing gap and heavy quasiparticles. The conditions for the appearance of the 
Andreev window are determined explicitly. The Andreev window leads to 
peaks in the conductance spectrum. The a =t peak should be clearly visible 
only near the upper critical field. This result can be contrasted with other 



non-self-consistent calculations existing in the literature [13|, |20jj. The cal- 
culated differential conductance as a function of external bias allows for an 
experimental verification of the results provided the FF state with an s-wave 
symmetry of the gap is detected. 
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Figure 6: (Colour online). Conductance plots for different Z values for (a) BCS state, (b) 
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the signals from both spin channels, and (d) for a 3D system at H a = 4.49 T w i/ C 2. Note 
that the distance between the peaks is twice the Zeeman energy 2h — 2g^BH a . 
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